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I.  Introduction 

When  a  steel  cylinder  is  filled  with  a  high  explosive  and  the 
explosive  is  detonated  the  fragment  weight  distribution  depends  on  the 
location  of  the  initiation  point,  or  points,  the  cylinder  diameter, 
the  wall  thickness,  the  type  and  treatment  of  the  steel,  and  the 
composition  and  density  of  the  explosive.  A  scheme  to  get  the  fragment 
weight  distribution,  analytically,  once  these  variables  are  specified, 
has  been  a  goal  of  warhead  research  for  over  thirty  years.  Two 
aspects  of  the  problem  are  considered  here,  namely,  the  construction 
of  a  one  parameter  family  of  fragment  weight  distributions,  and  the 
assignment  of  values  of  this  parameter  to  various  high  explosives,  for 
one  particular  steel  type,  cylinder  size,  and  initiation  point 
location. 

Two  sets  of  fragment  weight  distribution  data  were  used.  The 

first  set  (the  NWL  data)  is  from  work  at  the  Naval  Weapons  Laboratory, 

Dahlgren,  Virginia1,  on  the  effect  of  various  steels  and  treatments 

on  fragment  size.  It  consists  of  data  sheets  covering  116  sawdust  pit 

firings  of  Composition  B  loaded,  uncapped,  steel  cylinders  ranging  in 

diameter  from  ^0  mm  to  8  in.,  with  explosive/metal  mass  ratios  (C/M) 

between  0.1  and  0.1!.  The  second  set  of  data  (the  NOL  data)  is  mainly 

2 

from  a  1953  NOL  report  by  Solem,  Shapiro,  and  Singleton  ,  where  the 
effect  of  explosive  composition  on  fragmentation  was  studied. 

1  2 

The  examination  of  the  NWL  and  NOL  data  showed  that  the  commonly 

•s 

used  Mott  distribution-^  gives  a  reasonably  good  fit  only  in  a  central 
portion  of  the  fragment  weight  range.  There  are  more  small  fragments 
and  fewer  large  fragments  than  the  Mott  formula  predicts.  This  diffi¬ 
culty  is  overcome  here  with  a  more  complete  formula  which  involves 
only  one  independent  parameter,  conserves  mass,  and  uses  the  Mott 
distribution  where  it  applies.  The  measured  weight  distribution  is 
somewhat  dependent  on  the  fragment  recovery  technique.  In  the  NWL 
and  NOL  experiments  considered  here  the  ejected  fragments  were 
collected  in  sawdust,  after  crossing  an  air  gap.  This  was  followed 
by  magnetic  separation. 

The  analysis  here  makes  use  of  the  envelope  of  the  fragment  weight 
distribution  curves.  The  envelope  will  be  seen  to  have  a  central  part 
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in  the  construction  of  the  weight  distribution  formula,  in  finding 
the  optimum  value  of  the  independent  parameter  for  particular 
applications,  and  in  comparing  natural  fragmentation  with  preformed 
or  controlled  fragments. 
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II.  The  Weight  Distribution  Formula 

The  following  notation  will  be  used.  Let  m  be  the  fragment 
weight,  in  grains  unless  otherwise  specified  (7000  grains  *  1  lb). 

For  a  particular  metal  casing,  let  W„,  be  the  total  casing  weight  in 
grains,  and  let  W(m)  and  N(m)  be,  respectively,  the  weight  in  grains 
and  the  cumulative  number  of  fragments  of  weight  greater  than  m.  Let 

W  =  w/wT 

N  =  N/Wt  (fragments/grain), 

and 

IT  =  7000  N  (fragments/lb.) 

The  symbols  N  and  N(m)  have  the  same  meaning.  Similarly, 

N  «  N(m),  W  -  W(m),  W  *  W(m),  N  «  N(m). 

Note  that  to  conserve  mass 

dW  =  m  dN  .  (2) 

From  this 


where  Nq  =  N(0),  the  total  number  of  fragments.  In  normalized  form 


Here  NQ  is  the  total  number  of  fragments  per  grain,  and  NQ  is  the 
total  number  of  fragments  per  lb. 


The  Mott  Formula  for  the  fragment  weight  distribution  can  be 
written 


N(m)  =  Nq  exp  -(m/y)1/2  1  , 

L  J 
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where  Nq  and  y  are  constants  for  a  particular  weight  distribution. 

The  value  of  the  constant  Nq  in  (4)  is  consistent  with  (3),  i.e., 

N  »  N(0),  the  total  number  of  fragments  per  grain.  The  graph  of 
o  1/2  1/2 

(4)  is  a  straight  line  in  the  m  '  ,  log  N  plane,  where  -y  '  log  e 

is  the  slope  and  log  NQ  Is  the  ordinate  intercept. 

If  the  Mott  formula  is  assumed  to  hold  for  the  entire  weight 
distribution  then  the  constants  Nq  and  y  are  related.,  through  the 
conservation  of  mass,  by 

No  *  l/(2y) .  (5) 

To  show  this,  note  from  (4)  that  m  =  y[ln(NQ/N) ]2 .  Inserting  this 
into  (3)  leads  to  (5).  The  Mott  distribution  is  thus  seen  to  depend 
on  only  one  independent  parameter. 

The  cylinder  fragmentation  data  examined  here  shows  that  in  most 
cases  the  Mott  distribution  (4)  holds  only  over  a  central  portion  of 
the  weight  range.  There  are  more  small  fragments  and  fewer  large 
fragments  than  the  Mott  formula  predicts.  To  account  for  this  type 
of  distribution  a  three  part  formula  was  constructed.  Account  is 
taken  of  the  total  weight  of  fragments  weighing  less  than  one  grain, 
but  the  number  and  weight  distribution  of  these  fragments  are  not 
considered.  The  fragments  weighing  more  than  one  grain  are  divided 
into  three  groups.  The  Mott  formula  (4)  is  used  for  the  central 
portion  of  the  fragment  weight  range,  the  Mott  formula  with  the 
exponent  1/2  replaced  by  a  value  between  0.5  and  1  Is  used  for  the 
larger  fragments,  and  a  power  law  is  used  for  the  smaller  fragments. 
The  weight  distribution  is  written 


I  : 

N  = 

mp  ,  m^  <  m  <  m^  , 

ml 

=  1 

(6) 

II  : 

N  = 

r\zl  exp  [-(m/yI]:)1/2]  , 

m^ 

<  in  <  m3 

(7) 

III  : 

N  = 

nm  exp  [-(m/yjn)0*]  > 

m  > 

m3 

(8) 

4 
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The  subscripts  I,  II,  and  III  refer  to  regions  and  the  subscripts 
1,  2,  and  3  refer  to  the  region  boundaries.  As  in  the  Mott  formula, 
the  distribution  is  assumed  to  be  completely  determined  by  a  single 
parameter.  The  parameter  used  in  the  NWL  work  is  used  here,  namely, 

H,  the  average  weight  of  fragments  weighing  more  than  one  grain,  i.e., 

H  «  W(l)/N(l)  =  W(l)/N(l)  .  (9) 

The  notation  =  N(r^),  Wj  =  W(m^)  will  be  used.  Thus,  since 
m^  *  1,  the  left  boundary  of  Region  I  (Equation  6), 

N1  =  N(l),  Wx  =  W(l) ,  Nx  =  N ( 1 ) ,  Wx  =  W(l)  , 

but  N2  refers  to  the  cumulative  number  of  fragments,  per  grain,  of 
weight  greater  than  m2  grains.  In  this  notation  Ff  »  W1/N1. 

Since  the  weight  distribution  is  assumed  to  be  entirely 
determined  by  the  single  parameter  FT,  the  quantities  N^,  n^,  n^jj, 

WII*  yIII*  p*  q>  m2*  and  m3*  aPPearinS  in  (8)  -  (8)  are  functions  of 

ff. 

The  weight  fraction  of  fragments  weighing  less  than  one  grain  was 
found  *7om  the  experimental  data  in  reference  1.  The  data  were  fitted 
to  the  form 


V?T  -  Wx  K 

wT  FT  +  K  ' 

with  K  a  constant.  Then  from  (1),  (9),  and  (10) 


(10) 


-  M/OJ+K)  , 

and 

Nx  =  1/ (M+K)  . 


(11) 

(12) 


Here,  and  N.,  are,  respectively,  the  weight  fraction  and  number 

of  fragments  per  grain,  of  weight  larger  than  one  grain.  Note  the 
distinction  between  W^,  the  weight  fraction  weighing  more  than  one 
grain  and  Wj,  the  weight  fraction  of  fragments  In  Region  I. 
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When  the  weight  distributions  N  vs  m  are  plotted  for  various 
values  of  the  parameter  Ff,  the  curves  are  seen  to  have  an  envelope, 
determined  by  the  condition 

(3N/3?T)  *  0  (envelope).  (13) 

m 

Prom  another  viewpoint,  this  envelope  is  the  locus  of  the  maxima  of 
the  N  vs  Ff  plotted  for  fixed  values  of  m.  Denote  the  envelope  by 
N  (m  )  and  note  that  m  ,  and  consequently  N  ,  is  a  function  of  M. 

c  C  6  c 

1  2 

An  examination  of  available  experimental  *  data  shows  that  the 
envelope  points  (m  .  N  )  lie  in  Region  II,  or  at  the  II-III  boundary, 
for  all  values  of  Ff  considered.  We  thus  have,  from  (7)  and  (13),  on 
the  envelope 


1  dnII 


nII  ~  me 


1/2  dyII 


-1/2 


”dFT 


=  0  (envelope).  (15) 


Now,  (7)  holds  on  the  envelope,  so  that 


1/2- 


nII  *  N«  exp  3  . 


le  II' 


(16) 


Differentiating  (16)  with  respect  to  M  and  combining  the  result 
with  (15)  leads  to 


dN 


‘II 


1/2 


N 


(17) 


A  satisfactory  fit  to  the  experimentally  determined  envelope  had, 
for  M  larger  than  *1  grains,  the  form 

(M  >  il  grains)  (18) 


®e  =  CM  »  Ne  =  L/rae 


where  C,  B.,,  and  L  are  constants.  From  (17)  and  (18), 


-B2 


VII  =  (me/M)  =  K  (C/11)  M  , 


(19) 
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and  from  (16),  (18),  and  (19), 


-B, 


nII  =  ^Ly/C^  ^  exp[2] 


(20) 


Now,  from  (7)  and  (20), 

N2  *  nI;[  exp  [-(m2/pl;!;)1/2]  #  (21) 

and 

N3  =  nII  exp  [“(“j/Wjj)1^2]  .  (22) 

Also,  using  (6)  and  (7)  at  the  I-II  boundary,  the  exponent  p  in 
(6)  is 

p  =  lr.  (N2/N^)/ln  ra2.  (23) 

To  find  njTj,  q4a^d  PjTj  in  (8),  note  first  that,  since  (8) 
holds  at  the  II-III  boundary. 


nIII  s  N3 


exp  [(m3/yIIT 


(2H) 


where  is  known,  from  the  experimental  data.  Nov;  stipulate  that  N 
and  (3N/3m)jj  are  continuous  across  the  II-III  boundary.  This  leads, 
using  (7)  and  (8),  to 


UIII  =  £2c*ra3 


q-1/2  lt  1/2-  1/c* 
PII 


] 


(25) 


We  can  now  get  q,  the  exponent  in  III, by  conserving  mass,  i.e.,  by 
satisfying  (2).  Denote  the  normalized  weights  of  the  fragments  in 
each  of  the  regions  by  Wj,  WXI,  and  WIXI.  Prom  (17), 


^1  +  ^11  +  WIII  =  — 
x  x  M+K 


(26) 


The  weights  and  calculated  with  (2)  are 


N 


W1  “  T~ TT/p7[1  -  (VN1)(1  +  1/P>]  > 


(27) 
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Vr'n  *  N2  [m2  +  2(jn2MII)1/2  +  211.^]  -  [ir^  +  2(m3^III^1/2  +  2Pm]  • 


(28) 


Also,  using  (2), 


o 


‘III 


(In 


n 


III)1/q 


IT 


dN 


(29) 


Equation  (29)  was  solved  for  q  by  the  method  of  false  position,  with 
the  intergral  evaluated  by  Simpson^  rule  at  each  stage.  The  quanity 
Win  ava*lable  from  (26)— (28),  while  and  r.jjj  are  known 

functions  of  q  from  (24)  and  (25). 

It  turned  out  to  be  feasible  to  take  the  envelope  of  (7)  as  the 
II-III  boundary,  i.e.,  m«  =  m  ,  N-  =  N  . 

J  6  J  fe 

To  sum  up,  the  fragment  weight  distribution  (6 )— ( 8 )  considered 
as  a  one  parameter  family,  can  be  calculated  once  the  weight  fraction 
less  than  one -grain,  che  envelope,  and  the  region  separation  lines  are 
obtained  from  the  experimental  data.  The  distributions  need  be 
calculated  only  once  for  each  value  of  FT.  To  the  extent  that  the  one 
parameter  assumption  is  valid,  the  entire  weight  distribution  of 
fragments  weighing  more  than  one  grain  is  determined,  for  a  naturally 
fragmenting  cylinder,  by  specifying  the  value  of  FT. 

A  computer  program,  in  BASIC,  to  calculate  the  weight  distribu¬ 
tion  with  the  scheme  outlined  in  this  section  is  listed  in  Table  I. 
Table  II  is  an  outline  of  the  computer  program.  It  also  contains  a 
list  of  the  parameters,  the  corresponding  notations  in  the  text,  and 
the  equations  used  to  calculate  them.  Values  of  the  parameters  for 
the  region  divisions  and  the  particular  constants  L,  C,  B2,  and  K 
used  here  are  listed,  for  various  values  of  FT,  in  Tables  III  -  VI. 
Tables  III  -  V  can  be  used  to  construct  the  weight  distributions 
directly,  from  Equations  (6)- (8),  without  the  computer  program. 

Table  VI  shows  the  weight  fraction  in  each  of  the  regions,  for 
various  values  of  M. 
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III.  The  Fit  to  the  Data 

1  2 

Both  the  NWL  and  the  NOL  data  were  used  to  make  a  fit  to  the  form 
in  Equations  (6)- (8).  Each  data  sheet  furnished  by  NWL  covered  one 
sawdust  pit  firing  and  contained  lists  of  the  number  and  total 
weight  of  fragments  in  the  weight  ranges  0-1,  1-2,  2-5,  5-10,  10-25, 
25-50,  50-75,  75-150,  and  over  150  grains.  The  NOL  results2  appear 
as  plots  of  W(ra)  vs  m  for  0.5  gram  (7.7  grain)  weight  intervals 
starting  at  0.5  grams.  Each  plotted  point  in  the  NOL  report  is  an 
average  of  data  from  three  different  firings.  The  fit  here  was 
made  entirely  by  plotting  the  data  in  convenient  ways  and  drawing  best 
lines  by  eye. 

The  weight  fraction  of  fragments  weighing  less  than  one  grain  was 
found  from  the  NWL  data1  alone  (the  smallest  weight-  ran^  in  thn  NOL 
data  is  0. 5-1.0  gram).  In  Figure  1  the  NWL  data  is  used  to  plot,  on 
reciprocal  ruling  paper,  the  percent  by  weight  of  fragments  weighing 
less  than  one  grain  vs  M.  Note,  in  the  figure,  that  these  fragments  are 
11?  of  the  total  weight  when  M  is  5  grains,  6%  when  M  is  10  grains  and 
2%  when  M  is  30  grains.  The  data  was  fitted  to  (10),  with  K  *  0.6,  so 
that  from  (12), 

N.  =  7000/ (M+0. 6).  (30) 

The  envelope  (18),  constructed  to  fit  both  the  NWL  and  the  NOL 
data,  is 

me  =  2050/Ne  =  0.236  I1,77  (FT  >  U  grains).  (31) 

A  fit  made  for  the  boundary  line  between  Regions  I  and  II  gave 

5/4 

n2  =  l+(2/9) (M-l)  .  (32) 

Recall,  from  Section  II,  that  tne  envelope  will  be  taken  as  the  boundary 
between  Regions  II  and  III,  i.e.. 


(33) 
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The  remaining  quantities  in  the  weight  distribution  formula  (6)-(8) 
can  now  be  calculated  by  the  method  outlined  in  Section  II.  We  get 

yII*  nII*  N2*  N3  and  p  froni  '19)-(23)>  and  q,  n^^  and  u-j-j-j-  by  the 
numerical  solution  of  (2H) ,  (25)  and  (29). 

The  region  boundaries  (32)  and  (33)  are  shown  in  Figure  2.  The 
exponent  q,  used  in  Region  III,  is  plotted  vs  M  in  Figure  3.  Figure 
shows  the  weight  distributions  for  several  values  of  M  calculated  from 
(6)-(8)  and  (30)-(33)-  The  heavy  straight  lines  correspond  to  Region 
II  where  the  Mott  distribution  with  the  exponent  0.5  is  used.  In 
Figure  5  the  fit  to  the  weight  distributions  is  plotted  in  a  differe  t 
way,  as  N(m)  vs  the  parameter  M  '.‘or  fixed  values  of  m.  The  curve  for 
m  =  1  is  just  (30) . 

The  NWL  data  were  plotted  in  the  form  shown  in  Figure  5  by  taking, 
for  each  m  corresponding  to  the  end  of  a  weight  interval,  one  M  point 
from  each  data  sheet.  In  Figure  6  the  data  is  shown  in  this  form 
together  with  the  fitted  curves  in  Figure  5.  The  maxima  of  some  of 
the  fitted  curves  Is  a  bit  higher  than  a  best  fit  to  the  data  would 
give.  This  came  about  from  the  imposed  requirement  that  a  single  fit 
be  made  to  both  the  hWL  and  the  MOL  data. 
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IV.  Effect  of  Explosive  Composition 

The  NOL  results  consist  of  the  data  for  15  different  explosives 
treated  in  NAVORD  Report  29332,  together  with  unpublished  data  for 
eight  additional  explosives.  These  additional  data,  obtained  in  1952- 
1953,  were  taken  from  NOL  Project  Notebook  *<352.  In  all  cases  the 
explosives  were  loaded  into  AISI  10&5  steel  tubes,  9  inches  long,  2.0 
inches  ID,  2.5  inches  0D,  with  a  hardness  about  Rockwell  100-B.  The 
experimental  arrangement  is  shown  in  Figure  7.  Since  the  densities  of 
the  various  explosives  are  different,  the  explosives  are  being  compared 
on  a  volume,  rather  than  weight,  basis.  The  cylinder  breakup  will  vary 
with  the  method  of  initiation,  the  wall  thickness,  the  cylinder 
diameter,  and  the  type  of  steel,  as  well  as  with  the  explosive 
composition.  Hence,  the  values  of  M  obtained  for  the  different 
explosives  apply  only  to  this  system. 

To  assign  a  value  of  M  to  each  explosive,  the  weight  distributions 
were  first  calculated  with  (6)-(8)  for  values  of  M  ranging  from  *<  to 
50  grains  at  intervals  of  0.5  grains.  Then,  for  each  explosive,  the 
weight  distributions  corresponding  to  the  different  values  of  Ff  were 
superimposed  ur.  the  data  points  and  the  value  of  M  which  fit  the  data 
best  was  selected. 

The  MOL  data  for  the  various  explosives,  together  with  the  calcu¬ 
lated  weight  distribution  corresponding  to  the  assigned  M  values,  are 
shown  in  Figure  8. 

Results  for  all  the  explosive  compositions  studied  are  summarized 
in  Tables  VTI  and  VIII.  Mote  here  that  the  smallest  value  of  M,  10.5, 
was  obtained  with  (90/10)  BTNEN/Wax.  Among  the  explosives  practical 
for  military  use  Composition  B  and  Composition  A3,  both  with  M  =  12, 
are  at  the  low  end  of  the  range.  The  fragmentation  from  these  is  sig¬ 
nificantly  finer  than  that  from  cast  TNT(fl  =  17.5)  or  pressed  TNT(M  = 
22).  Fragmentation  from  some  of  the  aluminized  compositions  is  In 
between  that  from  Composition  B  and  TNT.  For  example,  M  is  13.5  for 
HBX-1,  15.5  for  H-6,  and  17  for  Tritonal.  Note  that  80/20  Ammonium 
Nitrate/TMT  produced  the  largest  average  fragment  weight,  with  H  *  *<5. 
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V.  Comparison  with  Preformed  Fragments 

The  concept  of  a  design  weight  is  useful  in  connection  with 
applications.  One  wishes  to  maximize  N(n.d),  the  number  of  fragments 
weighing  more  than  m^  grains.  For  each  design  weight  there  will  be 
a  value  of  M,  corresponding  to  a  weight  distribution,  for  which  N(m^) 
is  a  maximum  (set  m  =  in  Figure  6).  The  locus  of  these  maxima  for 
different  val  ;es  of  md  is  the  envelope  (13),  also  shown  in  Figure  6. 

The  relation  between  preformed  and  natural  fragmentation  for  a 
given  design  weight  is  shown  in  Figure  ?.  This  contains  a  -eg  N  vs  log  ir. 
plot  of  the  weight  distributions  from  natural  fragmentation  for  differ¬ 
ent  values  of  M,  calculated  with  (6)-(8).  The  envelope,  for  T\  bigger 
than  ^  grains,  is  a  straight  line  parallel  to  7000/m^,  the  nunf  er  of 
fragments  available  from  preformed  fragmentation.  A  comparison  of  the 
envelope  and  the  preformed  fragment  line  shows  that  the  upper  bound 
for  N(md)  from  natural  fragmentation  is  29%  (100x2050/7000  -  see  Sec¬ 
tions  II  and  III)  of  N ( ir.^ )  from  uniform  fragments  of  weight  m^.  This 
upper  bound  for  natural  fragmentation  is  solely  a  property  of  the  way 
the  steel  casings  break  up.  In  any  particular  Instance,  it  may  not  be 
attainable,  because  of  cost,  strength  or  material  availability  problems. 
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FIG.  3  THE  EXPONENT  q,  USED  IN  REGT 
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FIG.  4  FRAGMENT  WEIGHT  DISTRIBUTIONS  FOR  SEVERAL  VALUES  OF  ft  . 
CALCULATED  WITH  EQUATIONS  (6)  -  (8) . 
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FIG.  6-1  THE  NWL  DATA  (SUPPLEMENT  TO  REF  1).  THE  FITTED  CURVES  WERE 
CALCULATED  WITH  EQUATIONS  (6)  -  (8). 


CALCULATED  WITH  EQUATIONS  (6)  -  (8).  (CONT) 


CALCULATED 
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FIG.  7  ARRANGEMENT  OF  CHARGE  AND  INITIATOR  FOR  FIRING  OF  NOl  CHARGES 
(NAVORD  REPORT  2933) 


O  NAVORD  REPORT  2933 
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FIG.  8  CUMULATIVE  FRAGMENT  WEIGHT  DISTRIBUTION  (2  IN.  ID  -  2.5  IN.  OD 
CYLINDER  FRAGMENTATION  TEST)(CONT.) 
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AMMONIUM  NITRATE  -  TNT  (30/20) 


FIG.  8  CUMULATIVE  FR\GMENT  WEIGHT  DISTRIBUTION  (2  IN.  ID  -  2.5  IN.  OD 
CYLINDER  FRAGMENTATION  TEST)  (CONT.) 
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CUMULATIVE  FRAGMENT  WEIGHT  DISTRIBUTION  (2  IN.  ID  -  2.5  IN.  OD 
CYLINDER  FRAGMENTATION  TEST)  (CONT.) 
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FIG.  9  CUMULATIVE  FRAGMENTS  PER  LB.  VS  DESIGN  WEIGHT 
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TABLE  I 

LIST-BASIC  COMPUTER  PROGRAM  FOR  FRAGMENT  WEIGHT  DISTRIBUTION  CALCULATION 


45  DIM  M< 30) 

50  MC l)  =  1 

52  MC  2>=  1*5 

5  4  M<  3)  =  2 

5  6  MC  4)  =  3 

58  MC  5)  =  4 

60  MC  6)  =  5 

62  MC  7)  =  7 

64  MC8) =  1 0 

66  M  C  9  >  =  20 

68  MC10>=30 

70  MC 1  I )  =  50 

72  MC 1 2»=  75 

74  MC  1  3)  =  1 00 

7  6  MCI  4)  =  1  50 

78  MC 1 5>=200 

80  MC  I  6>=240 

82  MC  t  7)  =  480 

9  5  REM  FSQ034  CALCULATES  NB  VS  MBAR 

9  6  F0R  K 3=  4  T0  50 

9  7  C-0«  2364 

9  8  B2=  1  •  77 

99  L=2050/7000 

100  M  1=H  3 

110  Bl  =  C/4 

1 17  LET  U2=  B1*M 1  * B2 

120  LET  K“ *  6 

140  LET  M2=  1  +  C  2/9)*CMl- 1  >  t  C  5/4) 

153  LET  M4=4*U2 
155  M  3=M  4 
160  N6=L/4/U2 

162  N3=N6 

164  N7=N3*EXPCCM3/U2>**.5> 

170  N2=N7*EXPC  - CM2/U2)**  •  5> 

190  LET  Nl  =  l/CMl+K) 

200  LET  P=L0GC M2/N 1 ) /L0GCM2) 

210  LET  B=l*l/P 

220  LEI  R1=N1*C 1- CN2/N 1 ) t B)/B 

230  LET  S1=N2*CM2*2*CM2*U2>  ».5*2*U2> 

240  LET  S2=N3*CM3*2*CM3*U2)».5*2*U2> 

250  LET  R2=S1-S2 

260  LET  R3=M1/CM1+K)-Rl-R2 

265  GO SUB  1000 

267  U3=  C  2*  0*  U2*  *  • 5*M3**CQ-«  5)>**C 1/0) 

269  N8=N3*£XPCCM3/U3>**0> 

2  72  PRINT"MBAR="M1 
2  75  PRI NT**K=  "X#  "N 1  B=  "L*  7000 
2  78  PRINT"ME="M4#  "NEB="N6*  7000 
280  PRINT"M2='*M2»  "M3=  "M3 

282  PRINT"SQRCM2)="SQRCM2)»  "S0RCM3>  =  "SQRCM3> 

28  5  PRINT"MU2="U2»"SQRCMU2)="U2t. 5 

28  7  PRINT"N2B=  "N2*  7000*  "N 33=  "M3*  7000 

288  PRINT  ”Q="0 

290  PRINT"MU3="U3 

291  PRINT"ESN2="N7»"ESN3=t,N8 

29  n  PRI N  T"W 1 / WT= "R l » "W2/ WT= "R2» "W3/WT= "R3 
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TABLE  I  (CONT.) 

294  PRINT  "C  WT-WC  1 )  >/WT=”K/<Ml>K> 

29  6  PRINT 

298  PRINT'*N"»  “SOR(>*)'S  ,,N8‘' 

300  FOR  K2=  1  TO  17 

305  M=M<K2> 

310  IF  «<M2  THEN  500 
320  IF  M<M3  THEN  600 
330  N=N8*EXP<-<M/U3)**0> 

340  GO  TO  700 

500  LET  N=Nl*Mtp 

510  GO  TO  700 

600  N=N7*EXP<-(M/U2>**«5> 

700  PRINT  M»  SORCMXN* 7000 

710  NEXT  K2 

720  NEXT  K3 

725  GO  TO  9999 

1000  A=«2 

1010  LET  F4=0 

1020  Q1=.S 

1030  Q2= 1*5 

1040  LET  R=01 

1050  G0SU81 500 

1060  LET  F1  =  F 

1070  LET  R=Q2 

1080  GO SUB  1500 

1090  LET  F2=F 

1100  LET  Q3=Q1  +  C02-01)*ABS<F1>/<A3SCF:>*ABS(F2>) 
1110  GOSUB  1500 

I  120  LET  F3=F 

II  SO  IF  ABS<F4-F3XlE-7  THEN  1300 
1154  I  FABS<  F3>  <  l  E-  5  THEN  1300 

1  160  LET  F4=F3 
1170  PRINT 

1180  LET  Q5=03-A*(  02-01) 

1190  LET  Q2=Q3*A*(32-ai) 

1200  LET  01  =  05 
1210  GS  TO  1040 
1  300  L£TG=03 

1305  PRINT 
1310  RETURN 

1  500  LET  U3=t2*R*U2».5*M3t<R-.Si)t<t/R> 

1510  LET  Y=EXP(<H3/U3)  »R) 

1512  LET  Y=Y*N3 
1 51 S  LET  H=N3/Y/100 
1 520  S=0 

1530  F3R  J=1  TO  100  STEP  1 
1540  F=<-LOG<J*H>> *( 1/R) 

1550  IF  J=100  THEN  1590 

1560  IF  A3SUNTCJ/?)-J/2)<.  1  THEN  1610 

1  570  LET  $=S«-4*F 

1580  GO  TO  1620 

1590  LET  S=  S*F 

1600  GO  TO  1620 

1610  LET  S=  S+  2* F 

1620  NEXT  J 

1630  LET  F=H*S*Y*U3/3-R3 
1640  RETURN 
9999  END 
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Computer  Program  Outline 


INPUT 

*  Specify  values  of  FT  and  m  for  which  jT(m)  are  to  be  calculated. 
Insert  constants  L,  C,  B2,  and  K. 

CALCULATION  PROCEDURE 


rPor  specified  values  of  M 

Calculate  listed  parameters 


Program  Name  Text  Name  Equation  in  Text 


Rl, 


U2 

UII 

19 

M2 

m2 

32 

M4 

me(=m3) 

18 

N6 

Ne(=N3) 

20 

N7 

nII 

21 

N2 

N2 

12 

N1 

»1 

12 

P 

P 

23 

R3 

WI»  WII»  WIII 

26-28 

Q 

q 

24-29 

U3 

ym 

25 

N8 

nm 

24 

Print  parameter  values 
-►For  specified  values  of  m 

Calculate  !T(m)  6,7*  or  8 

Print  m,  m1/2,  N(m) 

■"••Next  m 


Next  M 


^TS3T 


(g*?3rejag&*ai 


4 

5 

6 
7 


8 
9 
JO 
1  t 
12 

13 

14 

15 

16 
1  7 
18 

19 

20 


22 

24 

26 

28 

30 

32 

34 

36 

38 

*0 


45 

50 

55 

60 

65 

70 

75 

80 

85 

90 

95 

100 


110 
120 
1  30 
1  40 
1  50 
160 
1  70 
180 
l^O 
200 
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REGION  BOUNDARY  AND  ENVELOPE  POINTS  FOR  VARIOUS  VALUES  OF  M 


ml 

m3 

N(m3) 

(GRAINS) 

(GRAINS) 

(GRAINS) 

(FRAGMENTS/GRAIN) 

1 

1*87738 

2*  74975 

*  106503 

1 

2*25708 

4.  03154 

7. 1 7516E-2 

1 

2*  661 5 

5*  63605 

S. 1961 4E-2 

1 

3*08673 

7. 40403 

3.95535E-2 

1 

3*  53023 

9*37814 

3*  1 2276E-2 

1 

3*98985 

1 1  *  552 

2.53512E-2 

1 

4. 4641 

1 3*9203 

2*  103R2E-2 

1 

4*95173 

16.4783 

1*  77723E-2 

1 

5*45173 

19*222 

1  *  52355E- 2 

1 

3*96323 

22* 1 476 

*013223 

1 

6*  48  551 

25*2519 

t*  159  74E-2 

1 

7*01  79  4 

23. 5313 

1 .02642E-2 

1 

7.  5599  7 

31  *9845 

9* 1562  IE- 3 

1 

3*11111 

35.  6075 

8  *  22458  E- 3 

1 

8*  67094 

39. 398  5 

7*  4332CE-3 

1 

9*23907 

43. 3552 

6.  7S433E-3 

1 

9*81515 

47. 4756 

6.  16« 53  E- 3 

1 

10*9899 

56* 1999 

5* 2 1 099  E- 3 

1 

12* 193 

65. 5573 

4.  467 19 E- 3 

1 

1 3*  4226 

75. 5354 

3*87709  £-3 

1 

14,677 

8  6.  1226 

3.  40047E- 3 

1 

1 5*955 

97. 3038 

3.00957E-3 

l 

17*2551 

1 09  *  03  4 

2*  68  468  E-3 

1 

18  *  S764 

121*  441 

2*41  152E-3 

1 

19.9179 

134*37 

2*  1  79  48 E-  3 

1 

21*2787 

1 47*865 

1  •  98058  E-  3 

1 

22*  658 

161*917 

1*808  6RE-3 

'  1828 

199*  45 

1.468  33E-3 

.09  3 

240*  339 

1  *213 52 E-3 

3  i29  7 

23  4*  505 

1  *  02936E- 3 

l 

">  3373 

331 *8  7S 

8*82431 E-4 

.  *2265 

382*  388 

7.  65S64E-  4 

45. 1925 

435*984 

6.7171 5E- 4 

49*231 1 

492*613 

5.9  449  7E-4 

53*  33.3  5 

552*226 

5*  3032 IE- 4 

57. 5115 

614.779 

4.  76362E-4 

61  *  747 

680  *  231 

4.  30526E- 4 

66*0425 

7  48  *  546 

3*9 1 23SE- 4 

70*  39  55 

319*  686 

3*  57230E- 4 

1  79.2655 

970*  31  5 

3*0181 7E- 4 

1  38*3416 

1131*37 

2*  58  736E-4 

1  97*6106 

1304. 15 

2*  24553  E- 4 

1  107*061 

1  43  6*9  4 

1  *  9  69  53E- 4 

1  116*683 

l  680*08 

1*  74312E-4 

1  126*463 

188  3*  39 

1 . 5549  5E- 4 

1  136*  403 

2096*  73 

1 . 39673E- 4 

1  1  46*  49  7 

2319*96 

1  *  26234E- 4 

1  156*727 

2552*94 

1*  1  4713E-4 

1  167*094 

2795*  57 

1*04753  E-4 

i 
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TABLE  ET 

CALCULATED  VALUES  OF  PARAMETERS  p,  M3r,  AND  N2 


M 

P 

nH 

(GRAINS) 

(GRAINS) 

(GRAINS) 

4 

-.531214 

•  78 69 S7 

.  63  7433 

- 

-•490195 

. 5301 76 

1.02039 

6 

- . 4541 43 

•  33  39  46 

1. 40901 

7 

-•437672 

•  292263 

1*85102 

8 

-. 429512 

. 230743 

2.  34454 

9 

-.425328 

•  187322 

2.838 

10 

-.423212 

. 1 55452 

3.  43  00  7 

1  1 

-•422237 

♦13132 

4. 1  1 9  53 

12 

-.421919 

•  1  12576 

4.B0S5 

13 

-♦421993 

9. 770S2E- 2 

5.53691 

1  4 

-.422302 

8. 569  41 E- 2 

6.31293 

15 

-.422754 

•  0758  43 

7.  1329  5 

16 

-.423292 

6.  76553 E-2 

7.99613 

17 

-.42388 

6.  0771 9E-  2 

3.90189 

18 

-  .  42449  4 

5  .  49  2  4  4E-2 

9.84963 

19 

-.4251 19 

4.991  18E-2 

10.8338 

20 

-•425745 

•04553 

1  1.8  689 

22 

- • 4269  78 

3«  8  S043E- 2 

14.05 

24 

-•428164 

3.  3003 3E- 2 

16.389  3 

26 

-.429293 

•  028648 

18.8333 

28 

-.430362 

2. 51 262E-2 

21. 5307 

30 

-.431373 

2. 2237SE-2 

24-3272 

32 

-.432329 

1.93373E-2 

27.271  1 

34 

-.433234 

1. 78183E-2 

30.3602 

36 

-.434092 

1 • 61 043E-2 

33.5925 

38 

-.434906 

J • 46346E- 2 

36.9662 

40 

-.435681 

1 . 33645E-2 

40-  479  4 

45 

-.437464 

1 . 08  49  5E-  2 

49  -  8  6  2  5 

50 

-. 4390S9 

9.00369 E-3 

60-03 43 

55 

-.440501 

7.60598E-3 

71-1262 

60 

-.441315 

6. 52033E- 3 

32.9639 

65 

-.44302 

5.  65901  E-3 

9  5.  59  7 

70 

-.4441 33 

4. 9  6334E- 3 

103.996 

75 

-.445166 

4.  39277E-  3 

123. 1 53 

80 

- . 446129 

3. 91857  E-3 

133  -  056 

85 

-.447031 

3.  S1936E-3 

153.69  5 

90 

-.4473  73 

3.  181  I8E-3 

170.053 

95 

-.448678 

2. 89  03  6E- 3 

137.  136 

100 

-.449  434 

2.  63996E-3 

204.921 

110 

-.450832 

2.2301 4E-3 

242. 579 

120 

- . 452102 

1.91  182E-3 

282.969 

130 

-•453264 

1 . 659 27 E-  3 

326.037 

140 

-.454335 

1. 45529E-3 

371.  ',36 

150 

-.455326 

•001238 

420.019 

160 

-.45625 

1.1  4S96E-3 

470.8  47 

170 

-.457113 

1 • 03205E- 3 

524.  182 

180 

-.457924 

9. 32749 E-  4 

579.989 

190 

-.45363  7 

8. 4762 4E- 4 

633.236 

200 

-.459  409 

7.  74059E- 4 

698.893 

N2 

(FRAGMENTS/GRAIN) 

. 1 S0747 
.  1  19813 
9.7,'  372E-2 
8* 0342  IE- 2 
6.  7 6 -*4 1  6E-2 
S.78262E-2 
S.00363E-2 
4.  38723E-2 
3*  88036E- 2 
•03461 1 
3. 1 I001E-2 
2.31P79 E-2 
2.SS371E-2 
2.33963E-2 
2. 1 4923E-2 
1.93261 E-2 
I.B3S8  5E-? 


1 • 5900SE-2 
1 . 39  326 E-2 
1 .2329 SE- 2 
1 . 10042E-2 
9-89  422E-3 
8  ■>  9  S420E-  3 
8.1501 6E- 3 
7. 45638 E- 3 
6.85301 E-3 
6*  32455E- 3 


5. 25660E- 3 
4.  45165E-3 
.003328 
3. 33373F- 3 
2.93454E-3 
2.  6069  IE- 3 
2. 33427E- 3 
2. 10466E-3 
1 .9092SE-3 
1 . 7  4 1  38  E-  3 
1  •  59  598  E-  3 
1 . 469  I  1  E-  3 


1 . 259 1 5E- 3 
1 • 09342E- 3 
9. 60039 E- 4 
8. 50908 E- 4 
7.60346E-4 
6.8  4267E- 4 
6. 19665E- 4 
5.6428  4E-4 
5.  16407E-4 
4. 74702E- 4 
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TABLE'S 

CALCULATED  VALUES  OF  PARAMETERS  n-jjj,  /i-jjj  ,  AND  q 


M 

nm 

(GRAINS) 

(GRAINS) 

4 

•  53801 5 

5 

•  338  579 

6 

*219969 

7 

*  1  549  41 

8 

•  1 1SS06 

9 

8*  9  7 608 E- 2 

10 

7*  19925E-2 

1  1 

5*9 1875E-2 

12 

4. 9 6358 E- 2 

13 

4* 2303  IE- 2 

14 

3*  65545E-2 

15 

3*  19476E-2 

16 

2*8 1969 E-2 

17 

2*  S0989E- 2 

18 

*022508 

19 

2*031 71E-2 

20 

1.84464E-2 

'“m  q 


(GRAINS) 

1*10108 

•  58  5232 

2*06463 

*644517 

3*32018 

*693007 

4*3394 

•73239 

6*60063 

*764517 

8*  58  769 

*790938 

10*  7882 

•812861 

13*1922 

*331207 

15.7917 

•846631 

18*  58 

*859327 

21*5511 

*871064 

24.  7 

*380722 

28*0221 

*839065 

31  *  51 33 

*8963 

35* 1699 

*902598 

38*9385 

*90309  7 

42*9  659 

•  9  1  29  1 

2<_ 

24 

26 

28 

30 

32 

34 

36 

38 

40 


1 • 54365E- 2  51  •  38  56 

1  •  31  38  4E-2  60.408 

1.13403E-2  70*0146 

9*90431 E- 3  80*1887 

8 • 73758E- 3  90*9155 

7*  77541 E- 3  102*181 

6*971 67E- 3  113*974 

6*  29266E- 3  126*282 

5.  71330E-3  139*095 

5*  21 457E-3  152*403 


*920341 
•926974 
*931727 
•935403 
*9  33229 
*940377 
*941976 
*943131 
*943922 
*94441  3 


45 

SO 

55 

60 

65 

70 

75 

80 

85 

90 

95 

100 


4.  23220E-  3  18  7*781 

3*  51  519  E-  3  226*073 

2-97422E-3  267*168 

2.SS49  4E-3  310*972 

2*  22271  E- 3  357*402 

1  *9  5448  E-  o  406*383 

1*  73447E-3  457*8  47 

1  *  551 53E- 3  51  1*  735 

1  •  39  7 58 E~  3  567*99 

1 • 26666E- 3  626*562 

1.15430E-3  68  7*  40 S 

1 • 05706E- 3  750*  475 


.944644 
*9  438  77 
.942466 
*94064 

*938  548 
*936289 
*933931 

•  931  521 

•  929093 

•  9266  7 
*924266 
*921894 


110 

120 

130 

140 

150 

160 

170 

l»o 

190 


8*  9  78  55E- 4  883*  135 

7*  7  3  783E-4  1024*26 

6*  75003E-4  1173*59 

5. 9  49  33E-4  1330*9 

5*  29  0  26E-  4  1  49  5*99 

4. 7  4052E-4  1668*67 

4.27666E-4  18  48*  78 

3*88 125E~  4  2036*  16 

3*  54114E-4  2230*67 


*917273 

*  9 1 28  4 

*908 609 

*90458 

.900746 

*897099 

.893627 

*890321 

*837169 
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TABLE'S! 


CALCULATED  WEIGHT  FRACTIONS  IN  VARIOUS  REGIONS 


M 

(GRAINS) 

l-W(l) 

WI 

Wn 

4 

. 130435 

9. 1 0696E- 2 

•  100636 

5 

.107143 

3*8  3233E- 2 

. 1 46453 

6 

9. 0909  1  £-2 

8 .90368  E- 2 

•176338 

7 

7.  39  47  4£-  2 

9.06! 19E-2 

•  197374 

8 

6.97674E-2 

9.22368 £-2 

•213024 

9 

•  0625 

v  36638 E-2 

•225162 

10 

5.  660 38 E-  2 

9. 48  3  62E- 2 

•  2  3  438  5 

1  1 

5.  1  7241 E- 2 

9 ♦ 57633E- 2 

•  2428  72 

12 

.047619 

9 • 64741 E-2 

•  249  5  7 

13 

4.  41  i  76E-2 

9 • 700 1 3E*  2 

•  255283 

14 

4.  109 59 E- 2 

9. 73752E-2 

•260223 

15 

3.8461  5E-  2 

•097622 

.264547 

16 

3-  61  446E-2 

9.  77638  E- 2 

•26837 

17 

3-  409  09  E- 2 

9 •  78 1 89  E-  2 

•  271779 

18 

3. 2258  IE- 2 

9-  78025E-2 

•  274843 

19 

3.06122E-2 

9. 77267E-2 

.277615 

20 

2.91262E-2 

.  09  7602 

•280137 

22 

2.  6543  7E-2 

9. 72377E-2 

•284566 

24 

2.  43902E-2 

9.6761  IE-2 

•288338 

26 

2. 2S564E- 2 

9. 62081 E-2 

.  29  1  598 

28 

. 0209  79 

9. 56033  E- 2 

•294451 

30 

1-9  60  78  E-2 

.09  4966 

.29  6973 

32 

1 .8  4049E-2 

9*  43076E-2 

.299223 

34 

•017341 

9* 36379 E-2 

•301246 

36 

1.63934E-2 

9*  29636E-2 

•303075 

38 

-015544 

9.22897E-2 

•  304741 

40 

1.4778  3E-2 

9« 1 619  7E-2 

•306265 

45 

1 • 31 579 E-2 

8.99775E-2 

.309  569 

SO 

1- IBS  #7E-2 

8-83987E-2 

•31231 

55 

1  •  079  1  4E-2 

8 • 68926E- 2 

•314629 

60 

9* 9 009 9 E- 3 

8 »  546 1 3E- 2 

•316621 

65 

9.  14634E-3 

8  -  41 029  E- 2 

•318356 

70 

8  •  49  8  58 E-  3 

8.28142E-2 

•319882 

7S 

7*93651  E-  3 

8. 15909  E-2 

•321238 

80 

7.  4441 7E- 3 

8 • 0428  6E- 2 

•3224S2 

85 

7.  00935E-  3 

7.9  3231 E- 2 

. 323547 

90 

6.  o22S2E-3 

7.82702E-2 

.32454 

95 

6*2761  SE-3 

7. 72661 E-2 

. 325446 

100 

5.9  642 1  E-  3 

7. 63073E-2 

.326277 

I  10 

5.  4249  5E-3 

7. 45127E-2 

•327748 

120 

4.9  751 2E- 3 

7.  28636E-2 

♦329013 

130 

4.  59418E-3 

7. 1 341 1 E- 2 

; 3301 16 

140 

4. 26743E- 3 

6.99  29  5E- 2 

•331086 

150 

3.93406E-3 

6.361  S4E-2 

•  331943 

160 

3. 73599 E- 3 

6.  73377E-2 

.  3 ->2  72 

170 

•  00351 7 

•066237 

•333416 

180 

3. 32226E- 3 

6. 51 552E-2 

. 334047 

190 

3.  1  479  5E-  3 

6. 4135SE-2 

•T 3 4622 

200 

2.99103E-3 

6.31717E-2 

•33515 

•  6778  6 

•  65808 1 
. 643716 
. 633067 
. 6249  72 
.618674 
.613675 
.60964] 
.606337 
• 603S99 
.601306 
.599  37 

.  S9  7722 
.5963!  1 

•  59  5097 
.  59 4046 
.  593135 


.  59 1  r  m 
.  590S1 1 
.  589637 
.  588966 
. 588  453 
. 5S8064 
.  58  7776 
.58  7  568 
. 5" 7425 
.58  7337 


. 587296 
. 587434 
.  58  7687 
.588017 
.  58S39  5 
. 588805 
.  589234 
.  589  6  7  5 
•  59012! 
.590  567 
. 591012 
.  $9 1  452 


.59  2  315 
.593148 
.  5939  49 
.59  4717 
.59  5453 
.596157 
.  59  683 
.  59  7  476 
.59809  4 
.  59B7-',8 


E 


[•  f> 


I- 

1 


i 


'& 


parameter  family  calculated  with  (6)-( 8) 
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